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The first 2 years of this program were used to make mid-IR observations of regions of star for- 
mation in the Orion nebula with the UCSD mid-IR camera at the UCSD/University of Minnesota 
telescope at Mt. Lemmon. These observations attempted to make the first systematic study of an 
extended region, known to have newly forming stars, and expected to have complex mid-IR emis- 
sion. We discovered, to our surprise, that most of the thermal emission originated from extended 
sources rather than from point sources. This interesting observation made the analysis of the data 
much more complex, since the chop/nod procedures used at these wavelengths produce a differ- 
ential measurement of the emission in one region compared to that in the adjacent region. Disen- 
tangling complex extended emission in such a situation is very difficult. 

In parallel with this work we were also observing comets in the thermal infrared, the other 
component of the original proposal. Some spectacular data on the comet Swift-Tuttle was 
acquired and published (Fomenkova, Jones et al., Astron. J., 1995, 110, 1866). This paper 
describes a changing jet structure observed over a 2 week period. The rotation period of the comet 
can be measured at 66 hours. The size of the nucleus can also be estimated (at 30 km) from the 
observed excess flux from the nucleus. These data have lead to the development of models 
describing the action of dust particles of differing sizes and composition leaving the nucleus. The 
spatial distribution of the predicted IR emission has been compared to the observed jet structures, 
leading to estimates of both particles sizes, relative amounts of silicate vs organic grains, and the 
amounts of dust emitted in the jets vs isotopic emission. This work forms the basis for the PhD 
thesis of James Sarmecanic. 
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ABSTRACT 

We report pre-perihelion observations of comet P/Swift-Tuttle obtained with the UCSD mid-infrared 
astronomical camera on 12 nights in November 1992. The images were taken through a 1 /in i wide filter 
centered at 11.65 /xm. During the observing run, the heliocentric distance decreased from 1.13 to 0.98 AU, 
and the geocentric distance increased from 1.16 to 1.33 AU. The spatial scale of the images was 700-800 
km per pixel. In addition, photometric data at wavelengths between 3.6 and 18.5 ^m were obtained on one 
night. The infrared images cover the cometary activity on time scales from hours to weeks and reveal large 
changes in the overall morphology of the coma. From periodic changes in the jet patterns w p e determined the 
period of nuclear rotation of 67.5 ±0.4 h. Photometry indicates that the temperature of the coma was 35% 
higher than the blackbody temperature at the same heliocentric distance and, hence, that the coma was 
dominated by small particles of average radius of 0.7 ^m. The dust mass loss rate varied with the 
heliocentric distance as /? -6 ' 1 * 3 and displayed —40% variation with the rotational phase. Two major jets were 
present in the images obtained on November 07-17, and a third area became active on November 24-29, 
increasing the average dust loss rate by a factor of 1.4. The relative positions of the three active areas on 
the surface of the nucleus are consistent with the positions of the most stable active zones identified by 
Sekanina [AJ, 86, 1741 (1981)] from the data on the P/Swift-Tuttle apparition in 1862. The radial brightness 
profiles suggest the radius of the nucleus is 15±3 km. This implies that the nucleus of P/Swift-Tuttle is —34 
times more massive than the nucleus of P/Halley. © 1995 American Astronomical Society. 


1. INTRODUCTION 

Studies of the nature and evolution of comets depend cru- 
cially on our understanding of the mechanisms and charac- 
teristics of cometary activity. During perihelion passage, 
comets are heated by solar radiation and release gas and dust 
into the inner solar system, thus forming a coma. From ob- 
servational studies, combined with extensive modeling, we 
know that the morphology of the coma is determined not 
only by the velocity and size distributions and the physical 
properties of the ejected material, but also by the geometry 
of jets and sources on the surface of the nucleus, its shape, 
and rotation geometry with respect to the Sun (review by 
Sekanina 1990). 

Temporal and spatial variations in the distribution of both 
gas and dust in coma indicate highly anisotropic emission 
mechanisms which are due to the inhomogeneity of nucleus 
surface properties and also depend upon the local insolation 
history. Gas molecules and very small dust grains display 
semichaotic motion (Horanyi & Mendis 1985), but the mo- 
tion of the optically efficient micron-sized and larger dust 
particles is much more collimated and predictable (Kitamura 


1987; Rabinowitz 1988). Due to the quasi-coherent behavior 
of dust particle jet flows, the trajectories can be dynamically 
extrapolated back to the surface of the nucleus, resulting in 
location and identification of their sources (Sekanina & Lar- 
son 1986). This makes studies of the dust features in the 
coma an indispensable tool for probing the nature of 
cometary nuclei. 

Imaging of comets in the mid infrared (3-20 /urn) is a 
very useful method for studying the morphology and the 
extent of the dust phenomena in the coma. The optical depth 
in the mid IR is very low (— 10 -5 ), which facilitates a 
straightforward interpretation of such images. The 10 /im 
wavelength range is dominated by thermal radiation from the 
dust itself and is sensitive to the dust distribution and its 
physical properties. A few 10 /xm images of comets have 
been acquired to date and have demonstrated the importance 
of this wavelength range for an insight into cometary pro- 
cesses (P/Giacobini-Zinner — Telesco et al. 1986; P/Hal ley — 
Campins et al. 1987; Hayward et al 1987; new comet 
Wilson — Campins et al 1989; new comet Austin — 
Fomenkova et al 1993). Analysis of these data has shown 
that the images are consistent with the premise that the over- 
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Table 1. Observational parameters for Comet P/Swift-Tuttle observations. 


Date (UT) 

Airmass 

Number 

of 

images 

R 

(AU) 

A 

(AU) 

Elongation 

(degrees) 

Phase Angle 
(degrees) 

Sun-Comet 
projection angle, 
North-to-East 
(degrees) 

Linear 

Resolution 

(km/pixel) 

1992 Nov 07.06-07.11 

L5 1-2.04 

12 

1.131 

1.166 

62.9 

51.1 

220.4 

702 

1992 Nov 08.07-08.11 

1.58-2.00 

8 

1.122 

1.165 

62.4 

51.3 

221.2 

701 

1992 Nov 09.09-09.11 

1.69-1.95 

4 

1.114 

1.167 

61.9 

515 

221.9 

703 

1992 Nov 12.06-12.11 

1.48-2.05 

9 

1.089 

1.174 

595 

51.7 

223.9 

707 

1992 Nov 13.07-13.11 

1.54-2.11 

12 

1.081 

1.179 

59.4 

51.7 

224.5 

710 

1991 Nov 14.06-14.11 

1.47-2.18 

14 

1.073 

1.184 

58.7 

51.7 

225.1 

713 

1992 Nov 15.06-15.07 

1.55-1-59 

2* 

1.066 

1.190 

58.1 

51.6 

225.7 

716 

15.08 

1.66-1.71 

2 b 







15.09-15.11 

1.90-225 

6 







1992 Nov 17.05-17.12 

1.44-251 

18 

1.052 

1205 

56.7 

51.5 

226.8 

725 

1992 Nov 24.07-24.11 

1.89-2.63 

8 

1.009 

1279 

515 

49.4 

2303 

770 

1992 Nov 26.07-26.11 

1.80-2.85 

10 

0.998 

1305 

49.9 

485 

2312 

786 

1992Nov28.05-28.il 

1.68-3.11 

15 

0.989 

1334 

483 

47.4 

232.1 

803 

1992 Nov 29.05-29.10 

1.67-2.77 

11 

0.985 

1349 

475 

46.9 

232.6 

812 


a. Taken with a 4.8 jj. filter. 

b. Taken with an 8.7 jx filter. 


all maximum activity takes place on the sunlit side of a ro- 
tating nucleus. In dusty comets with jets, the ejection of dust 
is restricted to a few discrete active areas representing only a 
small fraction of the nuclear surface. These areas become 
dormant when they rotate out of direct sunlight. A temporal 
lag in sublimation due to thermal inertia also should be taken 
into account when interpreting the diurnal evolution of the 
coma. 

In this paper we present 11-12 /xm images of comet 
P/Swift-Tuttle obtained pre-perihelion with the UCSD 
mid-IR astronomical camera. Comet P/Swift-Tuttle is a very 
active comet with a period of about 130 yr. It is the parent 
comet of the prominent Perseid meteor stream (Marsden 
1973). Apparitions of this comet in the inner solar system 
have been traced back to the years AD 188 and 68 BC (Yau 
et aL 1994). 

The observing run spanned the time period from 1992 
November 7 (IAUC #5654) to November 28 with data taken 
on 12 nights. We show an extended sequence of high- 
resolution images revealing the distinctive outlines of promi- 
nent features in the morphology of the coma. The temporal 
evolution of these features is diagnostic of the nuclear rota- 
tion. We also discuss how this rotation affects the activity of 
the comet. For one of these nights, we have photometric data 
obtained with the University of Minnesota (UM) infrared 
photometer at the O’Brien 0.76 m telescope. These data are 
used to derive the temperature of the dust, to estimate the 
characteristic size of dust grains in the coma, and for com- 
parison with other comets. 

2. OBSERVATIONS 
2.7 The Instrumentation 

Infrared images of comet P/Swift-Tuttle were obtained 
with the UCSD mid-IR camera (the “Golden Gopher”) on 
the UCSD/UM 1.5 m telescope at the Mt. Lemmon Observ- 
ing Facility. The telescope is optimized for observing in the 
infrared. At 1 1 fxm, the image spatial resolution is diffraction 
limited to X/D — 1.5 arcsec FWHM. The Golden Gopher 


platescale of 0.83 arcsec/pixel samples the point spread func- 
tion (PSF) with 2X2 pixels of the array. We used 64 
columnsX16 rows of the Si:As IBC array, which spanned 
53 X 13 arcsec. On warm nights or when the humidity is high, 
only a subset of the array can be used since the pixel satu- 
ration time due to higher background becomes less than the 
electronics readout time of 2.7 ms for a full frame. The cam- 
era provides good sensitivity in the wavelength range from 5 
to 27 /xm. At 11 /im, the noise equivalent flux density is 40 
mJy/arcsec 2 in 1 min (including both on and off source time) 
for a 1 /zm bandwidth. A detailed description of the instru- 
ment and observing procedure can be found in Pina et aL 
(1993). 

Photometry was carried out with the 15 filter multiaper- 
ture infrared GaGe bolometer on the 0.76 m IR telescope of 
the UM O’Brien Observatory. The bandpasses, calibration, 
and operational characteristics of the photometer are de- 
scribed by Ney (1974) and Hanner et aL (1990). Calibrations 
for the silicate filters were obtained by interpolating the cali- 
bration scales given by Gehrz et aL (1974, 1992). 


2.2 The Infrared Images and Photometric Data 

In November 1992, comet P/Swift-Tuttle was imaged at a 
wavelength of 11-12 )um during the 12 nights that weather 
permitted observations. On each night, from 4 to 18 images 
of the comet were acquired during 1-2 h of observing, the 
total number of images being 124. Table 1 summarizes the 
observing parameters and the Sun-comet-Earth geometry. 
During the observing run, the heliocentric distance R varied 
from 1.13 to 0.98 AU, and the geocentric distance A varied 
from 1.17 to 1.35 AU, resulting in a linear resolution at the 
comet of about 700-800 km/pixel (Table 1) or 840-960 km/ 
arcsec. A typical chop throw between the sky and the object 
images was 1 aremin, which was sufficiently large to avoid 
observable image contamination from a coma flux contribu- 
tion in the reference beam. 

If the dust in the coma is moving with velocity v d , then a 
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Fig. 2. The 2 to 20 /xm energy distribution of P/Swift-Tuttle. A weak silicate 
feature is superimposed on a smooth continuum caused by thermal emission 
of small grains. The fitted blackbody curve yields the grain temperature of 
360 K, which is 35% higher than the blackbody temperature at the helio- 
centric distance of 1.09 AU. 

crossing time required for the dust particles to cross the 
minimal resolution element is: 

T— A/ V d y 

where (p is an angular resolution and A is the geocentric 
distance. The crossing time defines an equivalent temporal 
resolution of the comet emission variability available from 
the images (Jewitt 1991). The terminal velocity u gas of the 
gas entraining dust grains into the coma can be estimated 
from the Bobrovnikov-Delsemme relation (Delsemme 
1982): 

D gas (km/s) = 0.58/r ,/2 (AU), (1) 

R being the heliocentric distance. Assuming that 
v d /v g z S —0.5-0.% (Finson & Probstein 1968) the temporal 
resolution is given by (Fomenkova et ai 1993): 

T(h) = 0.70<p(arcsec)A(AU)/?' /2 (AU). (2) 

Typical values for comet Swift-Tuttle observations are: 
<p=2X0.83" (two-pixel minimal resolution element), R— 1 
AU, and A=1.2 AU. The resulting temporal resolution 1 
h is comparable with the total observing time for each night, 
but is smaller than characteristic outburst duration of 2-3 h 
[as estimated by Sekanina (1981) for the previous apparition 
of P/Swift-Tuttle]. The typical time interval f 0 between suc- 
cessive images was 5-10 min, i.e., t 0 <T. Thus, images ob- 
tained on the same night were combined together to increase 
the signal-to-noise ratio and to extend the spatial coverage 
into the distant parts of the coma (Fig. 1, Plate 101). 

The photometric measurements of the IR energy distribu- 
tion of comet P/Swift-Tuttle were made on 1992 Nov 11.9- 
12.0 UT through narrowband filters at various wavelengths 
between 3.6 and 18.5 fx m (Fig. 2), including a number of 
data points in the silicate emission band at 8-13 /xm. For 


these measurements, the beam size and chop throw between 
the source and reference beams were 9 and 45 arcsec, respec- 
tively. Corrections, as described by Gehrz & Ney (1992), 
were made for a coma flux contribution to the reference 
beam. 


2.3 Flux Calibrations 

For absolute flux calibration, we used the stars a Aql, 0 
Peg, and a Tau. Calibrator images were obtained on each 
night, usually right before and right after the cometary ob- 
servations. The absolute stellar fluxes were taken from the 
IRAS Point Source Catalog (PSC, 1988) with a quoted un- 
certainty range of 4%-8% for flux measurement at each 
wavelength. However, since the satellite data are not color 
corrected and the IRAS bandpasses are much wider than our 
filter bandpass of 1 /x m, a correction factor needs to be ap- 
plied. This factor depends on the spectral index of the star in 
question. The flux densities, /(\ 0 ), quoted in the PSC for 
\ 0 =12, 25, 60, and 100 /xm, were derived assuming an in- 
trinsic spectral index of —1: 

or 

F=f q {\ 0 ) | /,(X)//,(X 0 )*x d\ 

=/,(X 0 )J (X/X 0 r'K* rfX, 

where F is the flux measured by a detector through a given 
filter, and R x is the relative filter response (PSC, Table 
II.C.5). If the actual spectral dependence of flux is then 

F=f(k 0 ) f /(X)//(X 0 )/? x d\ 

=/(Xo)J (X/Xor^x d\, 

and the actual color-corrected value of the flux density is 
/(X o)=/,(X 0 )/K, 

where the correction factor K equals 

K= J d\ j I (\/\ 0 y'R x d\ . (3) 

This correction has a large effect (up to 50%) and yields the 
true 12 /x m flux of the star. 

For our calibration stars, the spectral index (3 was esti- 
mated from the IRAS flux values at 12 and 25 /x m, assuming 
an average error of 6% for each flux measurement: 

/(12)//(25) = ( 12/25)“A 

IRAS spectra of those stars were checked in the IRAS Atlas of 
Low-Resolution Spectra (Olnon et al 1986) to make sure 
that they are featureless. To obtain true calibration values for 
our filter centered at 11.65 /xm, we applied an additional 
small correction (with the same index 0) to the color- 
corrected flux densities /( 12). The resulting fluxes for the 


i 


i 
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calibration stars are 24.2±1.5 Jy for a Aql, 286. 1 ±23.2 Jy 
for P Peg, 1 and 503.4±31.2 Jy for a Tau. 

Overall errors in the flux calibration of the cometary im- 
ages are due to (i) the intrinsic uncertainties in the IRAS 
measurements, as given above, (ii) sky noise in the calibra- 
tion star images, and (iii) airmass correction. Fluctuations in 
the sky emission are essentially random, and, since the raw 
sky flux is so large, the sky noise is very well described by a 
Gaussian distribution. Therefore, the uncertainty in the total 
flux in an aperture containing N pixels is equal to the average 
sky noise per pixel multiplied by the square root of N. On 
different nights, the sky noise contributes between 5% and 
12% to the uncertainty in the raw flux measurement for im- 
ages of a Aql, and about l%-3% to the uncertainty in the 
fluxes of p Peg and a Tau. These errors, root-sum-squared 
with the errors in quoted IRAS fluxes given above, yield 
uncertainties in the flux calibrations between 9% and 15% 
for a Aql, 7%— 11% for p Peg, and 5%-8% for or Tau. 

To account for the atmospheric extinction, we used an 
exponential decay approximation of the flux dependence on 
the airmass. For each night, the calibration coefficients ob- 
tained were plotted vs airmass on a log-linear scale. The 
slope of a linear least-squares fit to these data points yields 
an extinction coefficient a, which, on different nights, ranged 
from 0.04 to 0.25 mag per airmass. The uncertainties in the 
best-fit parameters of the airmass correction increase the er- 
rors in the flux calibrations. (Usually, either the number of 
calibration stars was too low or the airmass range covered by 
them was too narrow to yield a tight fit.) The final uncertain- 
ties in the flux calibration of cometary images are between 
9% and 18%. 

3. RESULTS 

3.1 Overall Morphology and Rotation of the Nucleus 

During its previous apparition in 1862, comet P/Swift- 
Tuttle was notted for its spectacular activity and was exten- 
sively observed [see Sekanina (1981) for the references]. The 
comet exhibited the dynamical evolution of three types of 
structures: jets, envelopes, and tail bands, indicating distinc- 
tive dust phenomena. Based on drawings and micrometric 
measurements of the jets made by the 19th century astrono- 
mers, Sekanina (1981) derived a nuclear rotation period of 
66.5 h and the obliquity of the orbit plane to the equator of 
80°, i.e., the rotation axis is almost in the orbit plane. He 
estimated that no more than 1% of the nucleus surface area 
was active, located eight discrete emission spots, and in- 
ferred a typical ejection duration of —2-3 h for each of the 
source regions. Based on orbital calculations, Yau et ai 
(1994) noted a surprising absence of nongravitational pertur- 
bations in the P/Swift-Tuttle motion over an interval of more 
than two millennia. They suggested that either the comet has 
a very large nucleus, more massive than typical cometary 
nuclei, or that its emission pattern has been symmetric and 
extremely stable during this period. We will show later that 
oui data support the notion of a very large nucleus. 


1 The error for (5 Peg is larger because this is an irregular variable with a flux 
variation of 5 % (Hanner 1994). 


Our infrared images cover the cometary activity on time 
scales from hours to weeks. The highly variable night-to- 
night activity of comet Swift-Tuttle is revealed in the coad- 
ded images (Fig. 1). Overall consistency of the images with 
an approximately 3 d period is immediately seen from the 
data and was noted by a number of other observers (Jorda 
et ai 1993; Feldman et al 1993). To estimate the rotation 
period more accurately, we sequenced the images vs rota- 
tional phase assuming periods between 50 and 72 h. The 
most consistent succession of jet development with time 
(Fig. 3) is reconstructed for a period T of 67<7<68 h. A 
period as long as 67 h, or as long as 68 h, changes the 
sequence in an unacceptable way. Therefore, our best esti- 
mate of the period is 7=67.5 ±0.4 h, which is slightly longer 
than the average period of 66.5 h derived by Sekanina (1981) 
from observations of the year 1862. 

The sequence of images on Fig. 3 reveals that a drastic 
change in the activity of P/Swift-Tuttle occurred between 
Nov 17 and 24. Two dust jets are distinguishable in the 
cometary images obtained between Nov 07 and 17. The first, 
strongest jet J1 is apparent in the beginning of the period 
(Nov 07): contours of constant brightness are elongated 
westward. It rotates counterclockwise (Nov 13), and dies 
away in the northwest direction by the middle of the rotation 
period (Nov 08, 14, 17). In these images the outer contours, 
determined by the pattern of earlier emission, are elongated 
to the northwest, but the more recent innermost contours are 
either circular, or elongated to the southwest. The images of 
Nov 14 and 17 show that the second, weaker jet J2 becomes 
active and is directed Sunward. This jet also rotates counter- 
clockwise (Nov 09), and disappears (Nov 12). At the end of 
the period (Nov 15) jet J1 is awakening again, pointing in the 
sunward direction. 

The pattern of emission from these two major active areas 
is in good agreement with observations of gaseous species 
CN, C 2 , and C 3 by Schulz et ai (1994) made on three nights 
in the beginning of October 1992. They estimated that, as- 
suming that the nucleus of P/Swift-Tuttle is in a state of 
simple rotation, the two active areas producing gas jets are 
separated by 130° on the surface of the nucleus. From the 
measurements of dust continuum at 0.483 and 0.523 pm 
Schulz et ai (1994) also suggested that the gas and dust jets 
in P/Swift-Tuttle are correlated in some way. Since our ob- 
serving run took place a month later than theirs, we cannot 
make a direct comparison of dust and gas emission. In our 
images, the time interval between the onset of jets J2 and J1 
is —22.5 h, which gives a distance between the active areas 
of 360° X (22.5 h/67.5 h)=120°. This result is consistent 
with the modeling by Sekanina (1981) based on the jet pat- 
terns observed during the previous apparition of P/Swift- 
Tuttle in 1862. He concluded that two most prominent emis- 
sion areas A and B were separated by 120° in longitude in 
cometocentric equatorial coordinates. The same areas may be 
active during the 1992 apparition and responsible for corre- 
lated gas and dust production! 

But the surface of P/Swift-Tuttle is not frozen perma- 
nently in the same, unchanging state. Images of Nov 24, 26, 
28, and 29 show the presence of a new, very strong, third jet 
J3. It is directed southward and has maximum activity in the 
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Fig. 3. Images vs rotational phase for a nuclear rotation period of 67.5 h. The brightest pixel of each image is marked by a cross. Contours are drawn at 60%, 
40%, 30%, 20%, and 10% level of the peak. The solid line through the brightest pixel of each image shows the sunward direction. 


beginning of the last quarter of the period (Nov 26), when it 
fully obscures the weaker jet J2. Then the activity of jet J3 
decreases 5 h later (Nov 29), while jet J1 reappears. The 
remnants of jet J3 are still seen in the beginning of the period 
(Nov 24) where outer contours of the coma are slightly elon- 
gated southward. The jet J3 is separated from jet J1 by —29 
h in time or by —155° and may correspond to the active zone 
C in the model by Sekanina (1981). In 1862, this zone be- 
came and stayed active closer to the perihelion than zones A 
and B. 

The sudden activation of a new strong jet indicates that 
the mechanical and morphological surface properties of 
P/Swift-Tuttle are heterogeneous and evolve with time. Se- 
kanina (1981) estimated that only about 1% of the surface 
area of the P/Swift-Tuttle nucleus was active during its 1862 
apparition. If the active areas responsible for the jets we 
observed are indeed the same which were active in 1862, 
then our observations give evidence of both diurnal and sea- 
sonal variations in the activity of the comet. The closer to the 
perihelion, i.e., to the cometary “summer,” the more areas 
on the surface change from a dormant to an active state. In 
Sekanina’s (1981) model of the 1862 activity of P/Swift- 
Tuttle, the maximum number of active areas occurs right 
after the perihelion passage, during the peak of the cometary 
summer. 

Alternatively, if the proximity of angular distances be- 
tween the three active areas observed in the 1992 and in 


1862 apparitions is a random coincidence, this would imply 
that the upper crust on the surface of the nucleus is unstable 
and subject to sudden openings of new crevices. Initiation of 
such a new active area would expose the inner, unmodified 
volatile-rich fresh material and result in a dramatic increase 
in the mass loss rate (see further). Fracturing of the surface 
may be caused by thermal stresses in cometary nuclei due to 
steep temperature gradients in the surface layer (Tauber 
& Kuhrt 1987). This explanation is also consistent with the 
observed jet activation at close heliocentric distance tf<l 
AU. 

3.2 Photometry and Comparison of Comet P/Swift-Tuttle 
with Other Comets 

Photometric data obtained on 1992 Nov. 12 were used to 
estimate the coma temperature (Fig. 2). A blackbody fit to 
the continuum emission (excluding the silicate band) yields a 
coma temperature 7 ob$ =360±10 K. The blackbody tempera- 
ture at a given heliocentric distance R can be determined 
from the energy balance between the power of absorbed so- 
lar radiation and the power emitted by grains in the thermal 
infrared: 

TTa 2 Q a (L Sm /4TrR 2 ) = 4T Tal Q e aT 4 0 ^ 

where Q a and Q e are the Planck mean absorption and emis- 
sion coefficients, respectively, a is the average grain radius, 
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33 Thermal Infrared Flux and the Mass Loss Rate 

In the steady-state model approximation of cometary dust 
emission, the dust production rate and the particle velocity 
do not depend on time, and particles are isotropically emitted 
from the nucleus. In this case, radial brightness decreases as 
Mr with the distance from the nucleus. The total thermal 
infrared flux F IR in a square aperture is roughly proportional 
to the mass loss rate M: 

Fi R =\n{\+V2){<p/A){3Q'(iT 4 ob J4Trpau ( i)M, 

or (6) 

F IR (W/m 2 )=1.16xlO' 25 [<p(arcsec)/A(AU)] 

x [ 7 1bs( K ) /c| (Mm)p(g/cm) 


Fig. 4. Comparison of P/Swift-Tuttle with other recent bright comets. The 
data on superheat and silicate excess in other comets are from Gehrz & Ney 
(1992). 

L Sun = 3.826X 10 33 erg/s is the solar luminosity, and 
a- 5.6696 X tO -5 erg/cm 2 /deg 4 is the Stefan-Boltzmann 
constant. Assuming Q a —Q e ~ 1, the blackbody temperature 
is 


7- bb (K) = (L Sun /167r C 7/? 2 ) 1/4 = 278.3/?- 1 ' 2 (AU). (4) 

Following Gehrz & Ney (1992), we define the superheat of 
cometary dust as 

S=T obs /T bb =(QJQ e ) m . (5) 

For the Nov 12 observations, R= 1.09 AU, T bb =267 K, 
5= 1.35 ±0.04, and Q e ^0.3<Q a . The high superheat is in- 
dicative of small grains present in the coma. If the grains 
were as large as, or larger than the wavelength of the emitted 
radiation 10 /xm, then their emissivity Q e would be close to 
unity, and their temperature would be close to that of a black- 
body [Eq. (4)]. The temperature has to rise above T bb in 
order for the grains to achieve radiative equilibrium. 

The emission in the silicate band (the four points between 
8.7 and 11.6 yum) exceeds the continuum fit by an average of 
0.34±0.15 mag. To compare the P/Swift-Tuttle dust compo- 
nent with other comets, we used the data by Gehrz & Ney 
(1992) (Fig. 4). The observed combination of the superheat 
and silicate excess characterizes P/Swift-Tuttle as an IR Type 
II comet. 2 Using calculations by Gilman (1974), as presented 
in Fig. 5 of Gehrz & Ney (1992), we estimated the average 
grain size a =0.7 ±0.1 yum. This is consistent with P/Swift- 
Tuttle classification as an IR Type II comet. For comet 
P/Halley, similar calculations confirmed that, usually, 
its coma was formed primarily by particles with radii of 
0.5-1 fim. 


2 The classification of comets into IR Type I and II was introduced by Gehrz 
et al (1989). Type I comets have continuum emission with a low superheat, 
weak or absent silicate features and Type I (ion) tails. Type II comets have 
10 and 20 yum silicate emission features superimposed on a superheated 
infrared continuum and prominent Type II (dust) tails. The observed differ- 
ences between dust emission can be explained by different grain sizes: large 
(5-10 ^m) particles in comets of IRType I and small (0.5-1 /zm) particles 
in comets of IR Type II. 


Xu^(km/s)]A/(kg/s), 


where (p is the size of a square aperture, a is the average 
grain size, p is the average density, and Q e =0.3. Obviously, 
comet Swift-Tuttle was not in a steady state as nonisotropic 
emission, temporal fluctuations in the dust production, and 
pronounced jets are clearly seen in the images. Many other 
parameters involved in these calculations are also very un- 
certain: the actual dust grain size, velocity and density dis- 
tributions, the dependence of Q e and T obs on the particle size 
and composition, etc. Most important is the dependence on 
the assumed particle size distribution. A considerable amount 
of the dust mass in the coma is in the form of large particles 
(>100 fim, McDonnel et al 1991), which we do not detect 
in 10 yum images (Campins et al 1989). Therefore, we can 
use the steady-state approximation only to obtain a rough 
measure of the lower limit of the dust mass loss. This esti- 
mate allows us a consistent study of short and long term 
variations in the dust activity of P/Swift-Tuttle during our 
observing run. Also, the dust production rate of other comets 
is usually estimated with similar assumptions, which makes 
the steady-state model a simple, but useful tool for compar- 
ing the productivity of different comets. 

The total thermal flux F m can be estimated from our mea- 
surements of the flux / u in a 1 yum wide filter centered on 
1 1.65 yum as 


IR' 




*12.15 


obs 


'11.15 


B bb( ^ > -^obs) td A 


(7) 


where A is the wavelength in /zm, and B bb (\,T obs ) is the 
Planck energy distribution. Assuming that on all dates the 
superheat of the cometary dust was the same as on Nov 12 
(i.e., the temperature of the grains was 35% higher than the 
blackbody temperature at the same heliocentric distance: 
7 0 b s =375.7 R~~ {n ) y we find that F IR =16.26 f n . The mass 
loss rate shown in Figs. 5(a) and 5(b) was calculated from (6) 
for a =0.7 yum, p=1.0 g/cm, v d -0A6R~ [f2 , and Q e =0. 3. 

Dust production in Fig. 5(a) varies approximately as 
F“ 63 , which agrees very well with the heliocentric depen- 
dence of the water production rate as derived by Hoban et al 
(1993) from a vast set of independent observations. This he- 
liocentric dependence is much steeper than a typical for com- 
ets R~ 2 -R - 4 one (Gehrz & Ney 1992). The increment of 
the mass loss rate at R<1 AU is caused by the activation of 
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a) Heliocentric dependence 



b) Rotational component 



Fig. 5. Dust mass loss rate of P/Swift-Tuttle. (a) Heliocentric dependence. 
The heliocentric distance was decreasing with time. The dotted line is the 
least squares approximation showing the /? -63 dependence of the dust mass 
loss rate, (b) Rotational phase dependence. The sinusoidal curve is used as 
the simplest periodic function. The width of each box is determined by the 
angular size of the corresponding image as follows from Eq. (2): the mass 
loss rate obtained from our measurements is an average over this time in- 
terval. The height of each box shows the error of measurement. 

jet J3 during the week between November 17 and 24 (cf. 
Sec. 3.1). Emission from the new area increased the average 
activity by a factor of 1.4. The overall activity is comparable 
with that of P/Halley 10 3 kg/s at R = 1 AU) and is about an 
order of magnitude higher than in other comets (Gehrz 
& Ney 1992). 

Day-to-day fluctuations in the mass loss rate [Fig. 5(b)] 
are roughly consistent with the suggested rotation period of 
67.5 h. On November 7-17, the dust production achieved 
maximum when the jet J1 was active, and dropped to the 
minimum during the intervals of the jet J 2 activity. The dif- 
ference between maximum and minimum dust production is 
~~40%, and similar variation in the water production was 
found from IUE observations (Feldman et ai 1993). The 
data presented by Schulz et ai (1994) show the following 
day-to-day abundance variations of gaseous species: 25% for 
CN, 26% for C 2 , and 40% for C 3 . This quantitative agree- 
ment confirms the hypothesis that, in contrast to comet 
P/Halley, the gas and dust jets in comet P/Swift-Tuttle are 
correlated. Remarkably, Feldman et ai observed a clear 
maximum in the water production on Nov 4. That is one 
rotation period earlier than our first observation of Nov 7, 
and, in accordance with the rotational dependence presented 
in Fig. 5(b), activity peaks at this rotational phase. 

3.4 The Size of the Nucleus 

Sudden activation of strong jets evidently contradicts the 
suggestion of a uniquely stable and symmetric pattern of 
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emission of P/Swift-Tuttle made by Yau et ai (1994). Se- 
kanina (1981) also inferred from his model calculations that 
different areas had been active during different rotations of 
the nucleus, the total number of active areas being up to 
eight. We consider the alternative hypothesis by Yau et al. 
(1994): that a surprising absence of nongravitational pertur- 
bations in the orbital motion of this comet is due to a very 
large nucleus. 

High-resolution infrared imaging gives us a key to direct 
observations of cometary nuclei even for highly' active, 
Halley-type comets near perihelion. Consider an ideal “unit” 
comet producing M = 10 3 kg/s of dust 3 in 1/zm size particles 
of the density p=l g/cm 3 at R = A=1 AU. In accordance 
with Eq. (1), the velocity of dust grains is v d ~0.5 km/s, and, 
with a superheat of 1.35, the temperature of grains is 376 K. 
From Eq. (6) we estimate that the infrared flux of such a 
comet observed in a 1 arcsec square aperture equals: 
/ =4.6X 10~ 12 W/m 2 . 

Suppose that the nucleus of the unit comet has an average 
radius Jt km (with a corresponding surface area of 4ttx 2 km 2 ) 
and radiates as a point source at a blackbody temperature 
T bb =278 K. The resulting infrared flux is given by 

C = <r7t b (*/A) 2 , 
or 

C(W/m 2 ) = 2.534X 10 _24 r£ b (K)jt 2 (km). (8) 

For je = 10 km, C equals to 1.5X10 -12 W/m 2 , implying that 
the contribution from the nucleus to the central brightest 
pixel is at least comparable with the flux from the dust. Note 
that the smaller the size of the central pixel (that is the higher 
the angular resolution), the larger is the relative contribution 
of the nucleus flux. 

In real images, however, the point source flux becomes 
spread out over some finite area described by the point 
spread function (PSF) of the instrument. Approximating our 
PSF with a Gaussian radial profile, PSF=(l/27r6 2 ) 
exp(— r 2 /2<5^), the contribution of the point source flux to the 
central pixel is reduced and given by 

Cpsp^ C/2 tt 6 2 , 

where S is determined by the full-width half maximum 
(FWHM) of the PSF as 

£ 2 = FWHM 2 /8 In 2 

for a Gaussian approximation of the PSF. Computer model- 
ing was applied to convolve the sum of the 1 fr dust contri- 
bution and the point source contribution from the nucleus 
with the instrumental PSF. The FWHM was estimated for 
each night from the calibration star images. The ratio of the 
nuclear flux C to the dust flux in the central pixel was used 
as a free parameter in the model. Values of C obtained from 
the best fit between the model and experimental radial 
brightness profiles yield an average nucleus radius of x = (15 
±3) km. This result is consistent with photometric measure- 


■bfhis value of M is close to the average mass loss rate in P/Halley at 1 AU 
and is about an order of magnitude higher than in most comets (Gehrz 
& Ney 1992). 
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ments by Williams who suggested that P/Swift-Tuttle 
nucleus should be at least 24 km across [reported in News 
Notes, Sky and Telescope 89(1), 1995, page 13]. For com- 
parison, the dimensions of the elongated comet P/Halley 
nucleus are (15.3X7.2X7.22) km, thus giving an average 
radius of 4.8 km (Szego 1991). 

We estimate that for the same density, the nucleus of 
P/Swift-Tuttle would be —34 times more massive than the 
nucleus of P/Halley. Hughes & McBride (1989) showed that 
the mass of a meteor stream is proportional to the mass of its 
parent comet multiplied by its dust-to-gas mass ratio. The 
mass of the meteor stream Perseid associated with P/Swift- 
Tuttle is (46±5)X10 12 kg and the average mass of the me- 
teor stream Orionid/Eta Aquarid associated with P/Halley is 
(1.5±0.4)X10 12 kg (Jenniskens 1994). Assuming the same 
dust-to-gas mass ratio in P/Halley and P/Swift-Tuttle, the 
ratio of masses of their meteor streams is 31, which is in 
excellent agreement with our result. 

4. SUMMARY 

Comet P/Swift-Tuttle at the heliocentric distance of 1 AU 
was one of the most active periodic comets after P/Halley 
and it demonstrated “Halley-type 1 ' behavior in some other 
respects too. Both P/Halley and P/Swift-Tuttle are old 
evolved comets having retrograde orbits (inclination >90°). 
Their comae are dominated by small (—1 /xm) particles. Im- 
ages of both comets were consistent with the presence of 
anisotropic flows of ejecta emanating from discrete active 
regions (vents) on the rotating cometary nucleus. Moreover, 
it was suggested that the same vents have remained active 
since their previous apparitions in 1910 for P/Halley 
(Rabinowitz 1988) and in 1862 for P/Swift-Tuttle (Schulz 
et al. 1994; this work). This hypothesis, if true, implies a 
certain degree of stability and, hence, might be used to infer 
the physical strength of the upper crust on the nuclei of old 
evolved periodic comets. Future missions landing on the sur- 
face of a cometary nucleus would shed more light on this 
problem. 

The activity of individual vents seems to remain stable 
from one apparition to another, but responds to the diurnal 
cycle and is affected by seasonal variations. To be fully ad- 
dressed in the future, it requires a long series of observations 
of a comet which would provide comprehensive information 
on day-to-day temporal evolution of cometary activity pre- 
and post-perihelion. Our observations of P/Swift-Tuttle, as 
well as' the interpretation of the 1862 observations by Se- 
kanina (1981), show that the number of active areas tends to 
increase closer to perihelion. In our case, there were two 
sources of jets on November 07-17, and the third vent be- 
came active in the images on November 24-29. The mass 
loss rate depended on the heliocentric distance as R ~ 63 . The 
initiation of strong emission from the new area increased the 
average mass loss rate by a factor of 1.4. From periodic 
changes in the jet pattern we determined the rotation period 
of P/Swift-Tuttle of 67.5±0.4 h. The dust production varied 
by —40% with rotational phase. 

In contrast to P/Halley, the dust and gas jets appear to be 
correlated in P/Swift-Tuttle. Observations of both comets 
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give rise to speculations about possible heterogeneity of their 
nuclei (Mumma et al. 1990; Fomenkova & Chang 1993; 
Hoban et al. 1993), but there is not enough data to draw 
more certain conclusions. 

The radial brightness profiles calculated for the mid- 
infrared images of P/Swift-Tuttle are consistent with the 
presence of a large nucleus, which contributes a non- 
negligible thermal infrared flux to the central part of each 
image. We estimated the radius of the nucleus to be 15 ±3 
km. If the densities of the two nuclei are equal, the nucleus 
of P/Swift-Tuttle is —34 times more massive than the 
nucleus of P/Halley. This may provide an explanation of the 
absence of nongravitational perturbations in the motion of 
comet P/Swift-Tuttle. 

We suggest that the differences in the general morphology 
of comae of new and old comets combined with the differ- 
ences in the dust grain size distributions may be the principal 
physical parameters correlated with the dynamical age of 
comets. For example, images of dynamically new comet 
Austin ( 1 989c 1 = 1 990V) obtained with the same UCSD 
camera did not reveal significant dust structures in the mor- 
phology of the coma, although the resolution w r as higher, 
<300 km per pixel (Fomenkova et al. 1993). Either the dust 
emission was more or less uniform from the w r hole nucleus 
(with a slight prevalence on the sunlit side), or the dust 
grains were ejected in a large number of microjets indistin- 
guishable as separate events from ground-based observa- 
tions. The long history of apparitions in the inner solar sys- 
tem causes formation of a thick dusty mantle, possibly glued 
together by complex nonvolatile organics, on most of the 
surface of a comet (Brin & Mendis 1979). Such an impen- 
etrable crust constraints the activity to a few active areas 
which are, probably, crevices reaching deeper volatile-rich 
layers. It effectively decreases and may be even completely 
prevents the dust emission from the rest of the surface. This 
scenario implies that the particle size distributions may be 
different in old and new comets. The coma of new comet 
Austin was dominated by larger particles (10 /xm) with the 
size distribution /(a) — 1 (Lisse 1992), while the coma of 

P/Halley contained a larger fraction of small particles (0.5-1 
jLtm) with the size distribution /(a)-fl“ 3 7 (McDonnell 
et al. 1991). The photometric measurements of P/Swift- 
Tuttle indicate that the average size of dust grains in its coma 
is —0.7 jam. A detailed modeling of dust emission processes 
in P/Swift-Tuttle to determine better its grain size distribu- 
tion is a subject of our future work. 
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